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Abstract 

Cavitation is a general phenomenon of the fluid flows with obstacles. It appears 
in the cooling conduits of the fast nuclear engines. A model of this phenomenon 
using the theory of Laplace and a common non-convex energy for the liquid and 
vapour bulks is proposed. This model makes it possible to determine a higher limit 
of the density of bubbles (a number of bubbles per unit of volume in the flow). The 
maximum intensity of cavitation is associated with the mechanical and thermal 
characteristics of the fluid flow. 
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1 Introduction 



The Laplace theory applied to a closed tank allows us to obtain an equilib- 
rium equation of a liquid in presence of vapour bubbles. We analytically prove 
that the total number of bubbles cannot exceed a maximum value. We extend 
the model to the mixture of fluid and gas in the case of steady flows without 
viscosity: the number of bubbles transported by the current admits an upper 
limit. 

It was observed in experiments that the number of bubbles created in a cavi- 
tating flow depends on the quality of the carrying fluid. It is well known that 
an injection of gas bubbles does not affect this upper limit. The number of 
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particles and microscopic gas bubbles is a significant cause of the intensity 
of cavitation [1-3]. Unfortunately, this parameter is generally inaccessible to 
measurement. 

The practical case is associated with the fiows existing in coolant circuits of 
fast nuclear engines; the used cooling agent consists of a mixture of sodium 
and argon. Sodium in liquid form or vapour is called "fluid" . Argon in form of 
gas, possibly dissolved, is a neutral constituent at a temperature much higher 
than its critical value. It is called "gas". It is not possible to analyze the qual- 
ity of the mixture. To estimate a maximum number of bubbles produced by 
a possible phenomenon of cavitation is very significant. It is an essential data 
of the nuclear security. 



2 Simplified model of cavitating flow 

We consider an isothermal permanent two-phase flow made up of fluid (for 
example sodium at the temperature of 400°C) and neutral gas far from its 
critical point (for example argon). We consider the approximation of a flow 
without viscosity. With an aim of simplifying the diagram of the system of 
cooling of fast nuclear engines, we separate the coolant circuit in two parts: 

- a first part of constant section s and volume v, in which the flow is car- 
ried along with cavitation. Cooling agent consists of sodium either in liquid 
form or in form of vapour bubbles; argon is dissolved in the liquid sodium 
or mixed with its vapour in the bubbles. The mixture contains n bubbles 
each one of volume vi . The bubbles are supposed in equilibrium in a reference 
frame convected by the flow and the edge effects of the pipe are not taken 
into account. The bubbles supposed without interactions are identical and in 
spherical form. 

- a second part of constant section sq and volume vo without cavitation in 
which sodium is in liquid form and saturated by dissolved argon. 

The motion of the mixture in each compartment is supposed to be uniform. 
One denotes respectively by u and Mq the velocity values in each compartment. 
The loop is closed and contains a total mass M of fluid (sodium). A pump 
discharges the fluid according to an imposed flow d. We denote by the 
average density of the fluid in the cavitating part. 

The theorem of Bernoulli is supposed to be applicable to each of the two com- 
ponents of the fluid-gas mixture. This simplifying assumption associated with 
a barotropic permanent motion makes it possible to connect the velocity to 
the density of each component. The theorem of Bernoulli extends to dynamics 
the equality of the chemical potentials between phases and components [4]. 
These relations can be replaced by more advancing equations of the motion 
representing, for a permanent flow, the connections between velocities and 
densities. 
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The volumic free energy of the fluid-gas mixture is supposed to be the sum of 
partial volumic free energies. The pressure of mixture is the sum of the partial 
pressures. The set of the assumptions leads to the following relations: 



Pm'" = Pi'" + ^^^1 {Pv - Pi) > 

which gives the value of the fluid mass in the zone of cavitation; pi and 
stand for the densities of liquid and vapour. It follows: 

(1) 

■"l^Pl-Pv^ 

The total fluid mass in the coolant circuit is: 

M = Pm" + PqVq, (2) 

where po denotes the density of the fluid in the part without cavitation. A 
pump gives the fluid flow d in the circuit, 

s p^u = sqPqUo = d. (3) 

Let us denote by F the volumic free energy of the fluid. The theorem of 
Bernoulli applied to each phase of the fluid component gives: 

^u^ + F' (pi) = ^ul + F' (po) (for the liquid phase) , (4) 
1 1 

-u^ + F' (p^) = -ul + F' (po) (for the vapour phase). (5) 

The approximations suppose a moderate cavitation. Although very synopsis, 
this schematization allows to get a first estimate of the maximum number of 
bubbles in the zone with cavitation. We take a state law of van der Waals type 
for expressing the fiuid partial pressure. The free volumic energy F is given 
by a single expression in the liquid and the vapour phases [5]. The relations 
(4) and (5) imply: 

F' iPi) = F' (pj (6) 
which is the equality of the chemical potentials of the two phases of the fluid 
component. The free energy of the fluid being a non-convex function of p, 
the values of pi and p^ are different. 

We denote by G the volumic free energy of the gas. We obtain G' (r^) = 
G' (ri) where is the density of argon in the bubble and r; the density of 
saturated gas dissolved in the liquid. Argon being a gas at a temperature far 
from its critical value, the function G is a convex function of the density r. 
Consequently, = ri. 

The partial pressures of the fluid and the gas are respectively: 

Ps{p)=pF'{p)-F{p) et Pa{r)=rG'{r)-G{r). 
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Within the framework of the Laplace model, the equilibrium of a bubble in 
the zone of cavitation is: 



(7) 



where a is the constant surface tension at the temperature of the flow and R 
is the radius of the bubble [6]. Owing to the fact that = ri and with K 
denoting (36%)^^^ a, the relation (7) is written in the equivalent form: 



F{p,)-F{p,)-F'{p,){p,-p,) = -Kv-, 



1/3 



(8) 



This relation remains unchanged if adding to F an unspecified linear function 
of p. Thus, F'{p) is chosen null for the densities of saturation pi^ and p^,g of 
the liquid and the vapour. 

The volumic free energy is supposed to be a regular function near the densities 
of saturation. The Taylor expansions to the second order in the vicinity of Pi^ 
and p^^ yield: 

FiPv)-\iPv-PvsfF"iPvs), 

F{Pi) = \{Pi-PisfF"{p,,). 

With taking into account the relation (6), we have 

F"{pJ{p.- pj^ F'\p,^){p,- Pi^), 
and the relation (8) yields: 



F"{Pis){pis- Pi){pi- Pv)^^ + \ 



1 Pis - Pi 



F"{p 



Is) 



Pi - Pv 



F"{pJ. 



1/3 



(9) 



The value C; of the sound velocity in the liquid is 



^ {pF' {p) - F {p))' 



Pis 



= PlsF" {Pls) . 



Pis 



The value of the sound velocity in the saturating vapour is 

Cl = PvgF" (p^g) . 

So, we obtain 

F" jPu) _ Cf p,, 

F"{pJ Qpu 
Far from the critical point of the fluid, we have 



F"{Pis) 
F"{pJ 



< 1. 
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Such is the case of hquid sodium at the temperature of the nuclear engine [7]. 
The range of values of and pi allows to write: 

Pis ~ Pi ^ 

Pi- Pv 

So, from the respective values of the various terms of relation (9), one deduces 
the approximated relation: 



3 An estimation of an upper bound value of the number of bubbles 



With the approximation (10), the relation (1) is written in the form: 

^ {Pl - Pm) 



n 



Pis 



By taking into account relation (10), one obtains the number of bubbles in 
the flow: 



n 



3Cf 



^ {Pls - Plf (Pl - Pm) ■ 



fill 



32na^Pi^ 

Let us notice that pi and are not independent but bounded by the relations 
(2), (3) and (4) expressing p^ as function of p^. Let us regard n as a continuous 
quantity; from these relations one deduces the differential system: 

V dp^ + vq dpQ = 0, 
u dp^ + p^ du = 0, 
uo dpQ + po duQ = 0, 

u du-\- F"{pi)dpi — Uq duo — F"{pQ)dpQ — 0. 
By eliminating du, duQ et dp^ between the four relations one obtains: 

(- + ^-F"{p,)^] dp^ = F"{p,)dp,. 

\Pm PqVo VqJ 

The respective orders of magnitude of the different terms of the preceding 
relation and the fact that Pq and Pi are close to Pi^ yield: 



-F"{p,)-dp^^F"{pi)dpi. 
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One obtains the derivative of as a function of Pi in the simphfied form: 



dpi V ■ 

Relation (11) where n is considered as a function only of p^, implies: 

l^'iPl) = 32^^3p^ ^ {~3 {Pls - Plf {Pl - Pm) + {Pls - Pit + ~)} ■ 

For Pl different of p^^, n'(p;) is null for: 

Pl-Pm = \{Pls-Pl)[^ + 'f)- (12) 

Let us notice that n = for pi = p,,. and Pi = p^. The value of Pi associated 
with relation (12) corresponds to a maximum of bubbles in the flow. The value 

of the number of bubbles per unit of volume is = — such as: 



Instead of the density of liquid, it is more realistic to consider the flow pressure 
in the part of the circuit subjected to cavitation. Indeed, this physical quantity 
is more accessible to measurements. 

Let us denote hj V = Ps{pi) +Pa{ry) and Vsat = Ps{pis) +Pa{rv): where Ps{Pis) 
is the saturated vapour pressure and Pairy) is the partial pressure of the gas. 
One obtains 

Vsat-V^Cf{pi,-pi). 

One deduces the value of the maximum number of bubbles per unit of volume 
in the form: 

N-^^^U^'^). (14) 



V 



For example, we consider the physical values of liquid sodium at the tempera- 
ture of 400 °C: surface tension o" = 0, 166 N/m ; sound velocity liquid sodium 
Ci = 2350 m/s ; density of saturating sodium pi^ — 856 Kg/m^ [7]. 

With choosing for the coefficient 1 H — ^ the value 10 and for difference of 

V 

pressures Vsat — V = 4000 Pa, one obtains in the part prone to cavitation a 
maximum of bubbles N — 1200 per liter . 
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4 Conclusion 

We presented an approach allowing to consider a higher limit number of the 
density of bubbles present in a closed loop of cavitating flow. This number is 
given by the relation (14). The effects of uncondensable gases (argon in our 
case) do not appear in the relation owing to the fact that the fluid (sodium) 
is saturated with gas. The bubbles are supposed in local equilibrium, which 
requires that their times of creation and collapse are weak with respect to the 
running time in the cavitating part of the circuit (a fine study could be carried 
out starting from the equation of Rayleigh-Plesset [8]). 

The relation binding the average densities of the fluid and of the liquid pi 
partially results from the relation (4). This relation corresponds to a perma- 
nent flow. It may be replaced by all other relations binding u,uo,pi et Pq. 
The method corresponds to a very simplified model of a flow of complex na- 
ture. Nevertheless, this model offers the advantage of showing that in a closed 
pipe, cavitation is a limited phenomenon. It should be noted that relation (14) 
takes into account the ratio of dimensions between the zones of cavitation and 
without cavitation. To decrease cavitation, one may reduce this ratio by creat- 
ing additional zones where cavitation is of no importance on the maintenance 
of the coolant circuit. 
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Abridged French version 
Resume 

La cavitation est un phenomene general des ecoulements fluides avec obstacles. EUe 
apparait dans les conduites de refroidissement des reacteurs nucleaires rapides. Un 
modele de ce phenomene utilisant la theorie de Laplace et une energie non convexe 

commune aux phases liquide et vapeur est propose. II pcrmet de determiner une 
borne superieure de la densite de bulles (nombrc do buUcs par unite de volume dans 
I'ecoulement). L'intensite maximum de la cavitation est associee aux caracteristiques 
mecaniques et thermiques de I'ecoulement fluide. 

Version frangaise abregee 

II a ete constate experimentalement que le nombre de bulles dans un ecoulement 
cavitant depend de la qualite du fluide. En particulier, le nombre de particules et 
de bulles microscopiques de gaz est un paramctre qu'il est important dc connaitre 
[1-2]. Malheureusement, ce nombre est difficilement mesurable dans la plupart des 
ecoulements. C'est notamment le cas pour I'ecoulement de sodium liquide en reacteur 
nucleaire de type RNR. D'autres experiences ont revele I'existence, pour un debit 
fixe, d'une valeur maximale pour le nombre de bulles; au dela de ce seuil, I'injection 
de microbulles de gaz ne semble plus modifier l'intensite de la cavitation [3]. 
L'objectif de ce travail est de quantifier l'intensite dc la cavitation a partir de 
I'cstimation analytique du nombre de bulles dans un ecoulement. Nous considerons 
un ecoulement diphasique permanent isotherme constitue d'un "fluide" (par exem- 
ple du sodium a la temperature de 400 °C) et d'un "gaz" neutre tres au-dessus 
de son point critique (par exemple de I'argon). Nous envisageons I'approximation 
d'un ecoulement sans viscosite. Dans le but de simplifler le schema du systeme 
de refroidissement des reacteurs nucleaires rapides, nous separons le circuit de re- 
froidissement en deux parties: 

- Une premiere partie de section constante s et de volume v, dans laquelle I'ecoulement 
s'effectue avec cavitation. Le refrigerant est constitue de sodium soit sous forme liq- 
uide soit sous forme dc bulles de vapeur; I'argon est dissous dans le sodium liquide 
ou melange a sa vapeur dans les bulles. Le melange contient n bulles, chacune 
de volume vi. Les bulles sont supposees en equilibre dans un repere entraine par 
I'ecoulement et les effets de bord de la conduite ne sont pas pris en compte. Les 
bulles supposees sans interactions entre elles sont alors identiques et d'une forme 
spherique. La zone de cavitation correspond a la region d'echange de chaleur entre 
les barres radioactives et le fluide refrigerant. 

- Une deuxieme partie supposee sans cavitation, de section constante sq et de vol- 
ume vo dans laquelle le sodium est sous la forme liquide et sature par de I'argon 
dissous. 

Le mouvement du melange dans chaque compartiment est pris uniforme. On note 
respectivement par u et uq les valeurs de la vitesse dans chaque compartiment. Le 
circuit est ferme et contient une masse totale donnee M de fluide. Une pompe en- 
traine le fluide suivant un debit impose d. Nous noterons la masse volumique 
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moyenne du fluide dans la partie cavitante. 

Le thcorcmc dc Bernoulli est suppose applicable a chacun des deux constituants 
du melange fluide-gaz. Cette liypotliese simplificatrice associee a un mouvement 
permanent et barotrope permet de relier la vitesse a la masse volumique de chaque 
constituant. Ces relations peuvent etre remplacees par des equations du mouvement 
plus cvoluees representant, lorsquc I'ecoulement est permanent, les liaisons entre la 
vitesse et la masse volumique de chaque constituant. 

Considerons pour simplifier des bulles en cquilibre dans un repere lie a cet ecoulement. 
En utilisant la theorie de Laplace, rencrgic libre totale du milieu s'ecrit comme la 
somme des energies de chaque phase et de I'energie superficielle des bulles. On ob- 
tient alors I'egalite des potentiels chimiques pour Ic constituant fluide et pour le gaz 
ainsi que I'expression du saut de pression a I'intcrface. Ces relations determinent 
les caracteristiques geometriques des bulles de vapeur. On represente I'energie libre 
du fluide supposee non convexe par deux paraboles dont les rayons de courbure aux 
sommets sont relies a la vitesse du son dans chaque phase. 

L'cnsemble des hypotheses efi^ectuccs conduit a I'expression (1) du nombre de bulles 
dans la zone de cavitation en fonction des densitcs du liquide et vapeur et de leur di- 
mension, a la conservation de la masse totale de fluide (relation (2)), a I'expression 
du debit impose (relation (3)). Le theoreme de Bernoulli applique a chacune des 
phases du fluide conduit aux relations (4) et (5) d'oii on deduit I'egalite (6) des 
potentiels chimiques [4]. De la definition de la pression a partir de I'energie libre et 
de la loi de Laplace on obticnt I'expression (8) reliant le rayon des bulles aux den- 
sites du liquide et de sa vapeur. En se limitant aux termes principaux, on montre 
que les deux relations (6) et (8) se reduisent alors a une seule expression: I'equation 
(10) reliant le volume des bulles a la densite du liquide dans la zone cavitante pi- 
La presence de gaz qui sature le fluide ne modifle pas la forme de I'equation de 
Laplace. 

Du fait des relations (2), (3) et (4), le nombre dc bulles n pcut etre considcre 
comme un fonction uniqucmcnt de p;. On constate alors que le nombre n de bulles 
admet une valeur maximale. Cette valeur prend une forme explicitable en fonction 
de I'ecart entre la densite du liquide et sa densite a la pression de saturation. C'est 
I'expression (13) donnant le nombre maximum de bulles par unite de volume. Elle 
constitue done une premiere maj oration analytique du nombre de bulles dans un 
ecoulement cavitant. 

La densite du sodium liquide n'etant pas une grandeur physiquement mesurable, 
il convient de transformer la relation (13). En utilisant une loi d'etat lineaire, on 

exprime le nombre maximum de bulles non plus en fonction de la densite du sodium 
liquide dans la zone de cavitation, mais en fonction de la pression totale. La rela- 
tion obtenue (14) montre que le nombre maximum de bulles par unite de volume 
est fonction de I'ecart entre la pression totale et la pression totale a la saturation. 
Un calcul numerique est effectue pour du sodium a la temperature des circuits de 
refroidissement des reacteurs nucleaires. II constitue un ordre de grandeur realiste 
du nombre de bulles des lors que Ton est capable de mesurer la pression dans la 
zone cavitante. 
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